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We report results of a search for particles with anomalously high ionization in events with a 
high transverse energy jet and large missing transverse energy in 2.4 fb _1 of integrated luminosity 
collected by the DO experiment at the Fermilab Tevatron pp collider. Production of such particles 
(quirks) is expected in scenarios with extra QCD-like SU(N) sectors, and this study is the first 
dedicated search for such signatures. We find no evidence of a signal and set a lower mass limit of 
107 GeV, 119 GeV and 133 GeV for the mass of a charged quirk with strong dynamics scale A in 
the range from 10 keV to 1 MeV and N — 2,3, and 5, respectively. 

PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 



One of the possible minimal extensions of the standard 
model (SM) is an addition of a new unbroken SU(N) 
gauge group [TJ [5] with new fermions, Q and Q, in the 
fundamental representation. Such a SM extension has 
two parameters: the new fermion mass vtlq and the 
strength of the new SU(N) gauge coupling, infracolor, 
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which can be defined through the scale A where it be- 
comes strong. 

If the new fermions, quirks, carry SM charges, they 
can be pair-produced at colliders. The production cross- 
section will scale with N, the number of colors in the new 
gauge group. In principle, the quirks could transform un- 
der any SU(N) group, but SU(2), SU(3) and SU(5) are 
considered explicitly. The phenomenology of such mod- 
els at hadron colliders in the case A <C toq ~ 0.1 — 1 
TeV has been recently studied in [3J, revealing an ar- 
ray of exotic final states that could have eluded previous 
searches for physics beyond the SM. The defining feature 
is that breaking of the infracolor string is exponentially 
suppressed due to the large value of the ratio tuq /A. Un- 
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like the SM quarks that immediately fragment into jets 
of hadrons, the quirk-antiquirk pair stays connected by 
the infracolor string like a rubber band that can stretch 
to macroscopic length 

L * 7? * 1 ^ ( iQO^ev) (lOOTev) ' (1) 

In this Letter, we consider a case where quirks have 
electric charge (e), no strong color charge, 10 keV < A < 
1 MeV, and 60 GeV < m Q < 1 TeV. This results in state 
with a decay length of many centimeters and a string size 
that is too small to be resolved in a detector, but large 
compared to atomic scales, e.g., for rriQ = 100 GeV the 
corresponding string size is less than 100 /xm, which is 
comparable to the resolution of the DO tracking system. 
The quirk-antiquirk pair will be reconstructed in the de- 
tector as a single highly ionizing track. Since the net 
charge of the string is zero, the track will be straight. 
Its experimental signature will resemble that of a single 
very energetic muon. However, due to low /3, the pair 
will often reach the muon chamber outside of the muon 
trigger time window. If it is produced in association with 
a high Et jet from initial state radiation, the event can 
be triggered and the track reconstructed. The signature 
which we consider is thus a high transverse momentum 
track with a large ionization-energy loss rate (dE/dx), a 
jet, and missing transverse energy (^?t) aligned with the 
track. 

Several collider experiments have searched for highly 
ionizing slowly moving particles, and set limits on pair 
production of stable sleptons [4], stop quarks [5] and 
charginos jU |7j- Although some of the searches [5] |S] 
were used to set an upper limit on the production cross 
section for a single, isolated, weakly interacting particle, 
the results can not be easily interpreted in the search for 
quirks. 

The present search is performed on data collected be- 
tween 2006 and 2008 with the DO detector [9] at the Fer- 
milab Tevatron pp Collider at ^/s = 1.96 TeV that cor- 
respond to 2.4 fb _1 of integrated luminosity. A detailed 
description of the DO detector can be found elsewhere 
[9TfTT]. It comprises a central tracking system in a 2 T 
superconducting solenoid, a liquid- argon/uranium sam- 
pling calorimeter, and a muon spectrometer. The track- 
ing system consists of a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT). The SMT, which ex- 
tends from radius of ~ 2 cm to ~ 10 cm, has a six-barrel 
longitudinal structure, each with a set of four layers ar- 
ranged axially around the beam and intersected with 16 
radial disks. The silicon wafers are 300 fim thick with 
strip pitches ~ 50 jim and ~ 150 fiia. In addition to 
position, the SMT detector provides a measurement of 
ionization energy with a dynamic range up to ~ 1.4 MeV 
per strip (~ 8 times the energy deposited by a minimum 
ionizing particle moving perpendicularly to the silicon 



wafer). The CFT, extending from a radius of ~ 20 cm 
to ~ 50 cm, has eight thin coaxial barrels, each support- 
ing two doublets of overlapping scintillating fibers. The 
muon system, extending from radius ~ 3 m to ~ 6 m, 
resides beyond the calorimeter and consists of a layer of 
tracking detectors and scintillating trigger counters be- 
fore 1.8 T iron toroidal magnets, followed by two similar 
layers after the toroids. The tracker readout gate is ca- 
pable of detecting slowly moving particles with velocity 
ft > 0.1 while in the muon system the readout gate is 
suitable only for particles with f3 > 0.3. The DO de- 
tector uses a three-level trigger system to select events 
for the offline analysis. The events used in this analysis 
were recorded using a jet trigger requiring in addition 
a substantial missing transverse energy calculated using 
the sum of the jet momenta (At = ISiets^H)- J e ^ s 
used in this analysis are reconstructed using the iterative 
midpoint cone algorithm [12] with a cone size of 0.5. Jets 
must satisfy quality criteria which suppress background 
from lcptons, photons, and detector noise effects. 

A modified [3] MadGraph [13] event generator with 
pythia [T3] fragmentation and hadronization is used to 
simulate signal events, which are further processed with 
a GEANT-based [T5] DO detector simulation and the same 
reconstruction as the data. We investigate quirk masses 
(tuq) in the range between 60 and 160 GeV and require 
that the jet from initial state radiation has a transverse 
momentum px > 65 GeV. In the simulation, quirks are 
treated as non-interacting particles and their energy loss 
in the SMT is calculated outside of GEANT. In the sce- 
nario considered in this study, the quirk-antiquirk pair 
forms a mesoscopic system. They follow a trajectory 
which is a classical superposition of the motion of the 
center of mass (c.o.m.) of the quirk-antiquirk system and 
the motion of the quirk (or the antiquirk) in the c.o.m. 
frame. The distribution of the velocities of the quirk- 
antiquirk system is very wide and peaks at (3 ~ 0.8 (0.2) 
for nig = 60 (160) GeV. In order to estimate the spe- 
cific energy loss dE/dx of the quirk-antiquirk pair, we 
integrate the energy deposition of each quirk along its 
trajectory and normalize the sum of both contributions 
by the corresponding displacement of the c.o.m. The 
superposition of both particles and such nontrivial kine- 
matics result in a large ionization (> 15 McV/cm) with 
a distribution that is substantially wider than that of a 
lepton or a hadron. To simulate the effect of detector 
resolution and instrumental noise, we smear the calcu- 
lated dE/dx with the width of the dE/dx distribution of 
muons from the Z boson decays measured in data. We 
simulate an effect of lost hits due to saturation that hap- 
pens at AE = J^l^i ~ IA Mcv wh z = 3Q is 

^ cos cos 8 5 f 

the thickness of a silicon wafer and where is the polar 
angle with respect to the proton beam direction and the 
quirk-antiquirk trajectory. 

The main SM backgrounds to the quirk signal are as- 
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sociated production of jets with a W boson and multi- 
jet events. We select candidate events that have ex- 
actly one jet with pr > 75 GeV and \rj\ < 1.6 [16] and 
|?t > 50 GeV. Events with additional jets with p^ > 
25 GeV are rejected to suppress multi-jet background. 
We require a track in the fiducial region of the central 
tracker with |^| < 1.6 andpx > 40 GeV. While the quirk- 
antiquirk pair is expected to follow a straight trajectory, 
we allow for a non-zero reconstructed track curvature and 
require pt > 40 GeV. The track should have at least four 
hits in the SMT and at least 12 hits in the CFT in or- 
der to suppress background from fake tracks. Multi-jet 
events are further rejected by applying the following iso- 
lation criteria. The sum of the px of all other tracks in 
a cone of radius TZ = J (A0) 2 + (Ar/) 2 < 0.5 around the 
candidate track must be less than 2.5 GeV. A similar 
isolation condition is applied for the total transverse en- 
ergy measured in the calorimeter cells in a hollow cone 
of radius 0.1 < TZ < 0.4 around the candidate track; this 
energy must be less than 2.5 GeV. We require the track 
to be back-to-back in the transverse plane with the jet, 
\A(f>(trk,jet)\ > 2.5. 

In order to suppress the W(— > tv)+jet background, we 
apply a cut on the azimuthal angle between the track and 
the missing transverse momentum, |A0 (trfc,$r)| < 0.5. 
We determine muon (electron) distributions of dE/dx 
from Z — > (e + e~) data events. In these events 

we require at least one of the leptons to have an associ- 
ated track that passes all signal selection criteria and an 
invariant mass of the two leptons consistent with the Z 
boson mass peak. The normalized distributions of dE/dx 
along such tracks are shown in Fig. [T] 

The background events where a candidate track is pro- 
duced by a charged hadron come from multi-jet events 
with incorrectly reconstructed $t- We construct the 
dE/dx distribution for these events from tracks found 
inside a jet (TZ < 0.5) with p T > 40 GeV and |r?| < 1.6. 
The track is required to pass the signal selection, includ- 
ing the track isolation requirement, and to originate from 
the same vertex as the jet (\Zj e t ~ z trk\ < 1 cm [H]). 

Another background comes from fake tracks recon- 
structed from hit patterns due to combinatoric ambigu- 
ities. To determine the shape of the dE/dx distribu- 
tion, we use the same event selection criteria as those 
used for hadrons, except the requirement of zero CFT 
hits to suppress contribution from real charged particles 
(see Fig. [I]). The shape of the dE/dx distribution for 
the fake tracks is significantly different from that for lep- 
tons and hadrons since crossing angles for such tracks are 
mismeasured and thus the thickness per layer crossed is 
incorrectly estimated. 

Fig. [2] shows the comparison between the data and 
the average of only the lepton and hadron backgrounds 
(without fake tracks), normalized to the region dE/dx < 
13 MeV/cm, and the predicted quirk signal for two mass 
values. We observe no excess of highly ionizing tracks 



above the expected SM contribution and set limits on 
the quirk production cross-section. 

We optimize the dE/dx cut to achieve the lowest 
95% C.L. expected limit for each quirk mass. The effi- 
ciency of the track quality requirement is obtained from 
Z — > fi + [i~ data events. Systematic uncertainties on the 
signal efficiency come from the combination of the un- 
certainties in modeling di-quirk dE/dx, trigger turn-on, 
and track quality, and are listed in Table I. We take the 
largest variation of individual lepton and hadron back- 
grounds from their average as a systematic uncertainty 
on the background. 

We observe an excess of data at low dE/dx with re- 
spect to the lepton and hadron-only background (see 
Fig [2]). A fit of the data to a sum of fake track and 
SM background improves the agreement at low dE/dx. 
However, in the absence of signal we disregard the po- 
tential contribution from fake tracks to the events being 
considered as quirk candidates. 

We perform a counting experiment in the region where 
dE/dx is larger than a ttiq dependent threshold, given in 
Table I, and set an upper limit on the quirk production 
cross section at the 95% C.L. Limits are calculated using 
the modified frequentist approach [17] and are consistent 
with expectation (Fig. [3] and Table I) . 

To summarize, we have performed a search for sin- 
gle highly ionizing tracks in events with an energetic jet 
and large missing transverse energy. This is the first 
study of this final state and the first search for quirk- 
antiquirk production. We find no excess of highly ion- 
izing tracks and exclude charged quirks of mass up to 
107 GeV, 119 GeV, and 133 GeV at the 95% C.L. with 
strong dynamics scale A in the range from 10 keV to 1 
MeV for SU(2), SU(3), and SU{5) gauge sectors, respec- 
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FIG. 1: (color online) Distribution of dE/dx for electrons, 
muons, hadrons, and tracks made of random hits (fake 
tracks). Each distribution is normalized to unity and the 
last bin contains overflow events. The shaded area represents 
both the statistical and the systematic uncertainties. 
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TABLE I: Results for each simulated quirk mass: dE/dx requirement, number of events in data, number of predicted background 
events, overall signal efficiency, quirk plus jet production rate, and observed and expected 95% C.L. upper limits on the signal 
cross section. 



M Q (GeV) 


dE/dx cut 
(MeV/cm) 


-/Vdata 


iVbkg ± syst 


Emciency±syst (%) 


N = 


^SU(N) 

2 N = 


(fb) 

3 N = 5 


Limit obs.(exp.) (fb) 


60 


19 


4 


5± 1 


11 ± 2.0 


63 


94 


157 


22 (25) 


80 


21 


2 


1.9 ±0.8 


9.9 ±2.1 


29 


43 


72 


20 (20) 


100 


24 





0.9 ±0.4 


9.2 ± 1.9 


16 


24 


60 


13 (18) 


120 


24 





0.9 ±0.4 


8.4 ±1.7 


9.3 


14 


23 


14 (19) 


140 


24 





0.9 ±0.4 


6.9 ± 1.4 


5.6 


8.4 


14 


18 (25) 


160 


24 





0.9 ± 0.4 


5.6 ± 1.1 


3.4 


5.1 


8.5 


22 (31) 
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FIG. 3: (color online) Observed and expected 95% C.L. lim- 
its on a{QQ+jet) for SU(2), SU{3) and SU(5) gauge sec- 
tors. The band shows ±1 standard deviation of the median 
expected limit. 
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FIG. 2: (color online) Distribution of dE/dx for candidate 
tracks in 4632 data events with overlaid background nor- 
malized to the number of events in data with dE/dx < 13 
MeV/cm. The shaded area represents the combined back- 
ground uncertainty. The quirk signal is shown for two masses, 
TfiQ — 80 (160) GeV, and normalized to the expected number 
of events, 13 (0.53). The last bin contains overflow events. 
The shaded area represents both the statistical and the sys- 
tematic uncertainties. 
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